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' A model for quantum dots is proposed, in which the motion of a few electrons in a three- 

. dimensional harmonic oscillator potential under the influence of a homogeneous magnetic field of 

arbitrary direction is studied. The spectrum and the wave functions are obtained by solving the 
classical problem. The ground state of the Fermi-system is obtained by minimizing the total energy 
' ^ ■ with regard to the confining frequencies. From this a dependence of the equilibrium shape of the 

(~| ' quantum dot on the electron number, the magnetic field parameters and the slab thickness is found. 
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PACS numbers: 73.20 Dx, 73.23.Ps 

I. INTRODUCTION 

Recent progress in semiconductor technology now allows the fabrication of artificially structured atoms in semi- 
conductors called quantum dots, in which electrons are trapped in a small localized region of space of a few hundred 



^ , Angstroms. In these nanostructures the electron wavelength is of the same length scale as the confinement so that 
' quantum effects are important, which results in the quantization of single-electron energy levels with excitation energy 

of several meV (see for review0~i) . A prominent feature of a Fermi-system which is confined to a finite region of space 

O . ... . Q 

is the emergence of shell effects in the single-particle spectrum as is well known in nuclear physicsu and more recently 



for large metallic clusters and references therein). In the present work we address this problem to mesoscopic 
systems like quantum dots which contain a small number of electrons. 

The electron states of few-electron quantum dots subjected to a strong magnetic field have been studied extensively 
The electrodynamic response of an interacting electron system in the presence of a confining potential is 
expected to be dominated by the many-body effects of the electrons. Sikorski and Merktl§ found experimentally the 
surprising result that the resonance frequencies in the magneto-optical spectrum are independent of the number of 
electrons in the quantum dot. The systems which are experimentally realized extent much less in the the z— direction 
than in the x — ?/— plane. This has led to a simple description of the far-infrared resonance (FIR) frequenciesEi~0 
as the energy levels of a two-dimensional harmonic oscillator potential in the presence of a magnetic fieldli. It was 
interpreted as a consequence of Kohn's theore nJi which has been generalized for a parabolic potentialE3~EIrE3~Erl. 
According to this theorem, the total Hamiltonian can be divided into two parts, the center-of-mass motion (CM) and 
the relative motion which contains the electron-electron interaction. Since the radiation of an external electric dipole 
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field couples only to the CM motion and does not affect the relative motion, the dipole resonance frequencies for 
the interacting system should be exactly the same as those of the nonintcracting system and be independent of the 
electron-electron interaction. 

A full understanding of the experimental results needs an analysis of many-body effects. Microscopic calculations 
using Hartree approximation for electron numbers N < loi neglected the exchange and correlation effects. The direct 
numerical diagonalization0 have been performed only for three- and four electron systems. The more complicated 
resonance structure observed intjtd raised the question as to the validity of Kohn's theorem for quantum dots. In 
order to describe the experimental data it was assumed that, in real samples, there is a deviation of the confining 
potential from the parabolic form, and different corrections have been introducedBEl. Also, the important role of 
the combined effect of the Coulomb forces and spin interaction leading to a redistribution of single-electron levels was 
demonstrated for the whole energy spectrum, especially for the low-lying single-electron states in the homogeneous 
magnetic fieldi!!!. The ground state transitions predicted in this way have been observed experimentallyEl. 

It seems therefore natural to assume that the properties of the electron states close to the Fermi level are determined 
by an effective mean field. It is true that the external field is the dominant part of the mean field, and thus the effective 
confining potential should reflect the main features of it. Yet it must also contain the effect of the interplay between 
Coulomb forces and the external fields which are governed by the charges in the adjacent layers and gates and the 
magnetic field. Due to these considerations, we assume that the confining potential should also take into account the 
changes that affect the properties of the single-electron states owing to a variation of the homogeneous magnetic field 
as well as the slab thickness. 

Based on the results discussed ir@ we conclude that the harmonic oscillator potential can serve as one of the 
phenomenological effective confining potentials in real samples, at least for small electron numbers. In this paper we 
study a simple model of a three-dimensional quantum dot in an arbitrarily oriented magnetic field for electron numbers 
for which the assumption of a parabolic type potential is still physically reasonable. Our procedure invokes an effective 
dependence of the parameters of the confining potential on the magnetic field. In this way, we accommmodate the 
different mechanisms mentioned above and we allow the system to adjust and change its shape under variation of 
the applied magnetic field and the particle number. We assume that the electron system behaves like an isolated yet 
confined droplet (see also discussion i One of the major result is a dependence of the single-electron states on the 
slab thickness. Such dependence shows in the FIR frequencies on which our analysis is focussed. We also report on 
transition strengths in absorption experiments and on the spatial extension of the confining potential in a quantum 
dot. Preliminary results have been reported in0. 
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II. THE MODEL 



We consider the Hamiltonian for TV noninteracting electrons moving in an effective mean field, i.e. a three- 
dimensional harmonic oscillator 

N N 

i=i j=i 

with A = [B X r\/2 being the vector potential associated with the homogeneous magnetic field B. The model does 
not explicitly take into account spin degrees of freedom. In other words, the system is assumed to be in a fixed 
spin state of an arbitrary degree of polarisation. Only the orbital motion is affected by the magnetic field. Using 
quantum mechanical equations of motion Yip0 calculated resonance frequencies for this model. Our approach is 
similar in spirit toi , but we first solve the classical problem, from which all quantum mechanical results easily follow. 
In contrast tc0'@ where the FIR frequencies are completely determined by the external potential, the effective mean 
field responds in our approach to the change of the magnetic field and particle number and thus affects the properties 
of single-particle excitations. We determine the frequencies uji, i = x,y, z hy minimizing the total energy for a given 
number of electrons N. In this way we allow the system to adjust the shape of the confining potential under the 
influence of the applied magnetic field and the particle number. 

The classical equations of motion for the cartesian components of the momentum and position coordinates read 

where p and r are combined to the six dimensional column vector denoted below as {p, r}. The matrix M is given by 

M = 
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where the notation fl = eB/{2mc) is used. Denoting the eigenvalues of A4 by ±iEk, k — 1,2,3 we obtain from the 
form A4 ~ UW the classical solution 

Mt), f(i)} = U cxp(I?i)VMO), f(0)} (4) 

where T> = diag(— — ii?2, —iE^, iE^, iE2, iEi). While the classical orbits are of little interest in the present context, 
the eigenmodes Ek and the system of eigenvectors listed in the matrix U are essential for the quantum mechanical 
treatment. 

The eigenmodes which are usually referred to as the normal modes are obtained from the secular equation 
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det\EI -M\=0 (5) 

with / denoting the six-dimensional unit matrix. The above equation turns out to be a third order polynomial in 
and has also been found bj0. The column vectors of 14 are the (complex) right hand eigenvectors of A4 (note that 
Ai is not symmetric, hence neither U nor V are unitary, yet V = U^^). We denote the column vectors of U by u^''\ 
they obey the equations 

{M + iEkI)u'-''^ ^ 0, 1,2,3 
{M - iEj^klW"^ ^ 0, A: = 4, 5, 6 (6) 

which can be solved as an inhomogeneous system by choosing for instance the sixth component of u'-'') equal to unity 
and then determining the normalisation as described in Appendix A. With the proper normalisation the classical 
Hamilton function of Eq. (|l|) is cast into the quantum mechanical form 

3 1 

H^Y.^E,iQiQ, + -). (7) 

The normal mode boson operators and Q are related to the quantised version of the classical coordinates by 

{p,r}=U{Q,Q^, (8) 

where we denote by a column vector which is the transpose of the vector (Qi, Q2, Qs, Q3, Q21 Qi)- The 

properties of the matrix W guarantee the commutator [Qk,Qk'] — iTT-Sk^k' as a consequence of [xi,pj\ = ih5i,j. We 
note that Eq. (j^) is the exact quantised equivalent of Eq. (|^) . 

III. QUANTUM ENERGY MINIMISATION 

The total energy of an A^-particle system associated with the Hamiltonian Eq.(0) is given by 

^^tot = 5Z^^'=K + l/2),- (9) 

The occupation numbers are the eigenvalues of Ql.Qk and take the values 0, 1,2,.. .. The ground state is determined 
by filling the single-particle levels J^k ^Ek{nk + 1/2) from the bottom. We take care of the electron spin only in 
obeying the Pauli principle which allows two particles in one level. It is clear that different sets of normal modes yield 
different sets of occupation numbers. The normal modes depend on the three components of the magnetic field and 
on the harmonic oscillator frequencies. From our assumption that the system adjusts itself under the influence of the 
magnetic field by minimising E'tot, a variation of the magnetic field strength leads to a corresponding change of the 
confining effective potential which is given by the oscillator frequencies. In other words, for a given magnetic field, we 
must seek the minimum of E'tot under variation of the oscillator frequencies. The variation cannot be unrestricted as 
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the confining potential encloses a fixed number of electrons, and assuming that the electron density does not change 
we are led to a fixed volume constraint which translates into the subsidiary condition uj^tOyiOz — with fixed. 
Denoting the Lagrange multiplier by A we solve the variational problem 

5{{g\H\g) - Xw^WyW,) = (10) 

where \g) denotes the ground state as described above. 

From Eq.(p^ we obtain, after differentiation with respect to the frequencies and using Feynman's theoremEl 

^{9\H\g) = {9\^\9), (11) 

the useful condition 

u;l{g\x'\g)^u;l{g\y'\g)^u;l{g\z^\g) (12) 

which must be obeyed at the minimum of -Btot- 

In this paper we restrict ourselves to consideration of a thin slab which extents essentially in two dimensions. This 
is achieved by varying only lu^ and ujy in the minimization procedure while keeping lu^ fixed at a value which is, say, 
five times larger than the other two frequencies. In this case only uj'^{g\x'^\g) = (^y{g\y'^\g) can be fulfilled. Choosing 
different (fixed) values of oj^ allows to study the dependence of the results on the slab thickness. 

It is known from the two-dimensional isotropic harmonic oscillator that shell closing occurs for Fermions at the 
numbers 2, 6, 12, 20, 30, .... At these numbers strong shell effects manifest themselves as the most stable quantum 
configurations. As a consequence, for Uz 3> lUx and i? = we must expect the minimum of Etot at the symmetric 
condition uj^ — ojy for such electron numbers. In the mean field approach, the breaking of spherical symmetry of 
the potential gives rise to a deformed shape (Jahn- Teller effeciS). For example, rotational spectra of nuclei and fine 
structure in the mass spectra of metallic clusters are explained as a consequence of deformed equilibrium potentials 
of these systems (se(£'i). Therefore, for numbers between the shell numbers a deformed configuration {uo^ ^ uoy) 
can be expected; we note that there are always two symmetric solutions: ujx > LUy and Wx < The condition 
ujz 3> LJx ensures that we have a genuine two-dimensional problem in that no particle occupies a quantum mode in 
the z-direction (n^ = 0). However, if the slab is made thicker {lOz smaller), the occupation of the first mode (n^ = 1) 
will occur at the top end of the occupied levels, which could mean that we find then the symmetric minimum at 
8, 14, 22, . . ., since two particles are occupying the z-mode that has become available. The distinction between a 
slab, that is sufficiently thin so as to prevent occupation of the first mode in the z-direction, and the slab which 
can accommodate the first mode is seen physically in the degeneracy of the two lowest normal modes at _B = 0. 
An estimate for lu'^ (see Appendix B) which is the frequency that just forbids occupation of a z-mode is given by 
> oj±{^/4:N + 1 — 3)/2 with uj± being the average of ojx and ujy. 
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IV. DISCUSSION OF RESULTS 



The level spacing Tiluq of the oscillator potential is determined by equating the Fermi energy ep of a free electron 
gas with the potential energy V{r) — ^m*LUQ < >. Assuming that the radius of a quantum dot grows with A^^/^ 
we have < 3/5i?QA^'^/'^. Consequently, we have chosen hcjo = 1.35A^~"'^/'^eF where the Fermi energy is obtained 

from the mean radius Ra and effective mass of typical quantum dots; for GaAs, i?o=320A and m* — 0.067me yield 
ep = 2meV. Throughout the paper we use meV as energy units, A for length and Tesla for the magnetic field strengths. 

A. The spectrum 

The three normal modes which can be discerned experimentally as excitation energies in FIR-spectroscopyQ behave 
in a distinctly characteristic way when the magnetic field is switched on. We recall that the values for uj^ and ujy are 
fixed by minimizing -Etot- In Fig- 1 we display typical patterns for two different electron numbers. In both cases, lu^ 
is chosen so large that remains zero. Note that, when the magnetic field is switched on, the occupation numbers 
ni,n2,n3 refer to the normal modes and can no longer necessarily be associated with an x-, y- or z-direction. 



E [meV] 




- B[Tesla] 

1 2 3 4 5 6 7 

FIG. 1. Normal mode frequencies as a function of magnetic field strength for A'^ = 16 (top) and A'^ = 30 (bottom). Solid, 
dashed and dotted lines illustrate results for 6 = O", lO" and 30", respectively. Values for huj^ are 9.2 in all cases. 

The degeneracy for A^ = 30 is clearly seen at _B = whereas three distinct modes, corresponding to a deformed 
shape, are seen at _B = for A^ = 16. The behaviour is shown for values oi 9 — O", 10° and 30". For 6 = the third 
mode (the highest at i? = 0) does not interact with the applied field, therefore there is a level crossing at about B = 5. 
This becomes a level repulsion ior 9 ^ as now all modes are affected by the magnetic field. As a consequence, the 
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^-dependence is strongly pronounced for the two upper modes, while such dependence is insignificant for the lowest 
mode. Here we find quite naturally an interpretation of avoided level crossing observed in the experiment^ associated 
with shape variations in a tilted magnetic field. We have also investigated the (/>-dependence of the normal modes, 
which is expected to become significant only for large values of 9] the differences are small, they can not be seen in 
Fig. 1. Below we return to a particular situation where the 0-dependence is important. 

E [meV] 
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0.5 1 1.5 2 2.5 3 
FIG. 2. Normal mode frequencies for TV = 30 and 6* = 0. The slab thicltness increases from the top to the bottom. 



There are windows of magnetic field strength for which the normal modes are virtually constant. They are discernible 
especially at smaller values of the magnetic field strength and are a reflection of shape changes. For = 30 the system 
is rotationally symmetric at i? = 0. At i? > 0.2 it changes into a deformed shape (w^ ^ LOy) which is associated 
with the larger derivatives of the two lower normal modes. The symmetric shape {uj^ = ojy) is now being restored in 
a continuous fashion when B is further increased, the two normal modes are virtually constant during this smooth 
transition. The constancy of the normal modes holds strictly only for 6 — Q and can be shown analytically; this is 
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deferred to Appendix C. At B = 0.7, where the system has eventually re-gained its symmetric shape, a second shape 
transition occurs for B > 0.7, which is reflected by the second sharp decrease/increase of the first/second normal 
mode. This pattern continues, yet it is less and less pronounced with increasing field strength. It remains to be seen 
whether these shape changes can be dissolved experimentally; finite temperature could disturb this pattern. 

Making the slab thicker, that is choosing Wz smaller so as to allow occupation of the next level of the third mode 
(n3 = 1), invokes distinct changes which could possibly be dissolved experimentally. In Fig. 2 we display the effect 
of varying the slab thickness for N ~ 30. For clarity we focus on smaller magnetic field strength and we display 
only results for d — '^i. The top figure is to facilitate comparison; the results are a repeat of Fig. (lb). In the middle 
figure the slab thickness has been decreased in such a way that the third mode becomes occupied when the field is 
switched on while it is unoccupied for zero field strength. As a consequence, the shape transition occurs for smaller 
field strength, it is here directly associated with the switch from na = to ns = 1 which is induced by the magnetic 
field in this particular case. The bottom figure refers to the case where = 1 throughout. The major effect is 
the lifting of degeneracy for B = 0. Since there are effectively only 28 particles occupying the two lower modes, the 
confining potential is deformed. 

With increasing field strength the potential changes into a symmetric shape which is attained at about \B\ = 0.55 
where it again undergoes a transition to a deformed shape in line with the discussion above. In this context we note 
that an increase of lOz beyond the values used in the results presented, which means an even thinner slab, pushes the 
upper level further away. A truly two-dimensional setting should therefore yield only two observable levels. 

E [meV] 




FIG. 3. Normal mode frequencies for N = 32 and ^ = 0. The slab thickness increases from the top to the bottom. 
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The effective occupation in the two lower modes can also be demonstrated by using iV = 22 or TV = 32. For 713 = 1 
this yields the effective numbers 20 or 30 in the planary motion with the expectation that the energy minimum is 
adopted for a symmetric shape. This is nicely confirmed and demonstrated in Fig. 3 where results for iV = 32 are 
presented. The top and bottom figures illustrate the cases where na = and 713 = 1, respectively. The interpretation 
is selfexplanatory after the previous discussions. 

In all the examples we see quantum mechanics nicely at work. Despite the simplicity of the model, these effects 
should be seen in an experiment. 

B. Dipole transitions 

The actual observation of the three modes is achieved by electromagnetic transitions, for instance in absorption 
experiments. Within the single-particle model considered here transition matrix elements between the ground state 
and single-particle excitations are available for any single-particle operator. The electric dipole transition is expected 
to feature prominently. We note that this transition rules out a spin fiip, which is in line with our assumption about 
fixed spin states. We interpret the three observed frequencieslll as the three normal modes. If no photon polarisation 
is measured in the initial and final states, the dipole transition strength is given by 

f,mf .m-y 

= E(l(/l^l5)l' + K/lj/l5>P + l(/l^l5>P), (13) 
/^™/ 

since the electric dipole operator is proprtional to f. Using Eq.(H) the matrix elements are given by a sum over the 
appropriate entries of the matrix U. The selection rules allow only transitions with ri^ ^ ± 1 by which all accessible 
final states are determined. Note that the configuration {n^} is different for different magnetic field strength. We list 
some results in Table 1. They refer to = and 713 = for the ground state. These figures do change as a function of 
9, electron number and slab thickness. The results demonstrate that there is a redistribution of dipole strength due 
to variation of these parameters and therefore of the shape of the system. This aspect could be the subject of further 
experimental analysis. 





S El 


S E2 


S E3 


B = 0.1 


1.0 2.0 


1.0 2.1 


0.7 9.2 


B = 1 


0.7 1.2 


1.2 3.5 


0.7 9.2 


B = 7 


0.07 0.3 


0.7 9.2 


1.1 12.4 



Table 1. Dipole strengths S for various values of the magnetic field. The associated frequencies are indicated next 
to the strengths. The numbers are normalised to the lowest frequency transitions at zero field. 
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C. Mean square values 

The calculation of matrix elements of and which are characterized by the shape of the effective confining 
potential follows similar lines as for the dipole matrix elements. In Table 2 we give a few values for = 30 and 
9 = Q. Note the strong deformation aX B — 0.8. For larger field strengths the mean values decrease, since the effective 
frequencies of the confining potential are augmented by the Larmor frequency; this effect also enhances a symmetric 
shape to an increasing extent for increasing field strength. We stress that these values all obey the condition Eq.(|2|). 
For 9 > Q the pattern remains essentially the same except for the fact, that the shape is not symmetric for large field 
strength but deformed; the degree of deformation depends on 9. This is again understood from the Larmor frequencies 
which are now different for the different directions. This argument is nicely confirmed when using 9 > and (j) — 45*^; 
now the symmetry is guaranteed by the choice Bx = By^ and, in fact, the system settles down to a symmetric 
configuration for large field strength. For B = 7,9 = 30", = the mean values are 1015 and 1069, respectively, 
while for </> = 45° they are both equal to 1044. Note that such pattern for large field strength is independent of the 
particle number. For instance, irrespective of the strongly deformed shape for = 16 at i? = 0, we find at B — 7 
also a symmetric shape for (p = 45*^ but a slight deformation for 9 = 30" and cj) = 0. Again we stress that the figures 
depend on the slab thickness in the spirit of the previous paragraph. 







V{9\y'\9) 


5 = 


1167 


1167 


S = 0.8 


770 


1773 


B = 7 


981 


981 



Table 2. Mean square values of the extension in the x — y— plane of the quantum dot under variation of the magnetic 
field strength. 

V. SUMMARY 

In view of the good qualitative agreement with experimental dataQ we think that the model gives a fair account 
about the essential features of quantum dots on semiconductor interfaces. The simplicity of the soluble model allows 
consideration of an arbitrarily oriented magnetic field. The minimal energy requirement invokes a re-arrangement of 
the quantum configuration which is manifested as shell effects and hence shape changes of the quantum dot. By this 
mechanism we find an effective A^-dependence of the FIR frequencies in quantum dots (see Fig. 1). An interesting 
aspect of the model is the constancy of the FIR frequencies (the normal modes) under variation of the magnetic field 
as long as the potential is deformed. As is discussed in Section 4 the system responds to the increasing magnetic 
field strength by changing the confining potential towards a spherical shape rather than by changing its energy. This 
is shown analytically in Appendix C for 9 = where it strictly holds. It is interesting to note that for 9^0 the 
statement still holds to high accuracy. As long as the shape of the effective confining potential is similar for different 
electron numbers, the splitting of the FIR frequencies due to the magnetic field should be equally alike. When the 
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tilting angle 9 is switched on, all three modes are globally affected by the magnetic field leading to an avoided level 
crossing or anticrossingO. A thorough experimental analysis of a variation of slab thickness, tilted angle and strength 
of the magnetic field, which leads to the manifestation of quantum effects such as shape variations for fixed and 
different numbers of electrons, could assess the reality and limitations of the model. 

The authors acknowledge financial support from the Foundation for Research Development of South Africa which 
was provided under the auspices of the Russian/South African Agreement on Science and Technology. R.G.N, is 
thankful for the warm hospitality which he experienced from the whole Department of Physics during his visit. 



Note added in proof: 

After submission of the present paper a report about recent experiments, which confirm essential aspects of this 
investigation, has appeared in Physical Review Letters, Vol. 77, 3613 (1996) by S. Tarucha, D.G. Austing, T.Honda, 
R.J. van der Hage and L.P. Kouwenhoven. 

APPENDIX A: NORMALISATION OF THE EIGENVECTORS 



The Hamilton function of Eq. can be written in matrix form 

H = {p,r}^H{p,r} 



(A.l) 



where 
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We aim at the quantum mechanical form 
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(A.4) 



Exploiting the fact that 
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(A.5) 
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where / is a 3 by 3 unit matrix, it follows that, up to normalisation factors, the matrix V = can be written as 
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(A.6) 



This implies that U^HU is in fact skew-diagonal as in Eq.(A.4), and therefore U can be normalised such that 
U'^nU = Uqn,. Using Eq.d), Eq.(^) follows. 



APPENDIX B: ESTIMATE OF uo' 







Consider the energy of the three-dimensional oscillator 

E = fiuJz{n^ + 1/2) + huj±{n± + 1) = E'^ + 



(B.l) 



where tOx 



UJy 



LUj_ ^ ujz- The condition uj^ S> 0Ji_ ensures that we have a genuine two dimensional problem 



in that no particle occupies a quantum mode in the z-direction [ux ^ Q,ny ^ 0,nz = 0), i.e. the first level with 
(rix = 0, Uy = 0, = 1) lies higher than the ones which are filled for given particle number. From the condition 

E^i^x = n-y = 0,nz — I) > E{nx + Uy = nj^ 7^ 0, = 0) it follows 



(B.2) 



For a two-dimensional system the number of particles (two particles per level) is yl = {N + 1)(N + 2) where N 



Uy = n±^ is the shell number of the last filled shell. Therefore, from the equations above it follows 



^ = (?ij_ + l)(rt_L -I- 2) i.e. n±_ 



%/4ATT-3 



Substituting Eq.(B.3) into Eq.(B.2), we obtain 



-,">^(x/4ZTT-3) 



(B.3) 



(B.4) 



APPENDIX C: CONSTANCY OF THE NORMAL MODES 



For 9 = Q the normal modes are 



El = ^(-^ 



El = ^(-^ 
E?. = uj,. 



-ml 
4r!2 



y^(.;2_^2)2+8J72(c.2+,.2) + 16f^4) 
y^(c.2_,.2)2+8J72(c.2+..2) + 16f^4) 



(C.l) 
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The last equation shows that the largest mode E3 is independent of the magnetic field strength, since is kept fixed. 
From the first two equations we obtain a relation between two sets of frequencies luix , ujiy and lu2x , '^2y which yield 
the same normal modes for different field strength. With the notation = fi^^ — > 0, where iliz and refer 
to the different field strengths Bi and B2, we obtain from the requirements Ei 2{Bi) = Ei^2{B2) 



1 



-2,-2 

To ensure that uj2x and uj2y are real, A must obey the condition 



4A2 + " ^ly? ~ ^^'Hx + ^ly) + 16A4) 



4A2 - ^{ojI^ - u:%? - 8A2(u;L + u;%) + IBA^) 



(C.2) 



(C.3) 



This means that the normal modes must change with the magnetic field strength for a spherical shape {loix = ^ly) 
as is established by the results. In turn, it is possible that Ei and E2 do not change under variation of B as long as 
LUix 7^ ^ly Further, if the magnetic field strength has increased up to the value where condition ( |C.3D becomes an 
equality, then the spherical shape {uj2x = ^2y) is attained. Recall that throughout this transition from deformed to 
spherical shape the normal modes (Eq.( p.l| )) have not changed. 

It remains to show that this solution is in fact the minimal energy solution. For 9 = the explicit expressions for 
{g\x'^\g) and {g\y^\g) read 



2m 



El E2 



El - El 



Si S2 
El E2 



(C.4) 



{9\y^\9) 
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2m 



Si 
El 



S2 
E2 



■417? /El 



Ef-El \Ei 



S2 
E2 



where = X], ("-fc + denotes the sum over all occupied single-particle levels. Suppose the condition Eq.(|T^) is 



fulfilled for the confining frequencies loix and uiy for the value of the magnetic field Bi{fliz) . Substituting Eqs.(C.4) 
into Eq.(|l^) and taking into account that ujix ^ ^ly, we obtain 

ml 



Si ^ Ti2 
El E2 



'ly 



El -El 



El E2 
El E2 



= 0. 



From Eqs.(C.2) it follows 



'^Ix + ^ly = + ^ly - 4(^^L - ^l,)- 



(C.5) 



(C.6) 



Using the result of Eqs.(C.4),(|C.5|) and (C.6), and the fact that Ei and E2 remain the same, we obtain 



u:lMx^\g2) - uj%{g2\v^\g2) = 
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where \g2) denotes the ground state referring to the frequencies L02xT^2y and the magnetic field B2. This means that 
Eq.(|l|) is fulfilled. 
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